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ABSTRACT 


The  effect  of  initial  spinal  alignment  on  the  location  and  mag¬ 
nitude  of  maximum  vertebral  stress  during  ejection  was  studied  using 
the  Orne-Liu  discrete  parameter  model  of  the  spine.  Face  curtain, 
shoulder  harness,  and  seat  back  restraints  were  added  to  the  model  as 
linear  springs.  Spinal  alignment  data  used  were  from  x-rays  of  5th, 
40th  and  95th  percentile  (sitting  height)  men  seated  in  the  MK-J5(D) 
ejection  seat  under  static  conditions.  Maximum  normal  stresses  were 
shown  to  occur  at  LI (5th) ,  T12(40th)  and  T9 (95th)  with  face  curtain 
and  shoulder  harness  restraint.  These  locations  correspond  almost 
exactly  to  the  predictions  of  injury  based  on  static  observations  of 
the  curvature  of  the  initial  configuration  of  the  pilot's  spine.  In¬ 
clusion  of  posterior-anterior  seat  back  support  decreased  maximum 
stresses  as  did  the  use  of  an  improved  lumbar  pad  which  placed  the 
lower  spinal  column  in  extension.  Failure  to  utilize  the  face  curtain 
restraint  gave  rise  to  large  normal  stresses  in  the  upper  thoracic 
column.  Results  indicated  that  a  state  of  nearly  uniform  axial  stress 
exists  in  the  column  during  ejection  and  thus  the  location  of  maximum 
bending  stress  dictates  the  spinal  location  of  maximum  normal  stress. 
Hence,  initial  spinal  alignment,  in  terms  of  the  curvature  of  the 
column,  is  a  major  determinant  of  the  location  and  magnitude  of  maxi¬ 
mum  normal  stress  for  a  given  set  of  restraints. 


ROBERT  W.  BAILEY 
Colonel,  MSC 
Commander 
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INTRODUCTION 


The  advent  of  high  speed  and  performance  aircraft  during  the 
early  1940 's  necessitated  the  development  of  powered  pilot  extraction 
systems.  The  ejection  systems  are  capable  of  separating  a  man  from  an 
aircraft  without  his  striking  the  vertical  stabilizer  and  attaining 
sufficient  height  to  allow  parachute  deployment  during  low  level  ejec¬ 
tion.  In  order  to  accomplish  this  during  the  stroke  length  of  the 
ejection  device,  the  man-seat  unit  must  be  subjected  to  a  relatively 
high  amplitude,  short  duration  acceleration  pulse  (Figure  1).  Unfor¬ 
tunately,  this  generates  large  forces  and  deformations  in  the  spinal 
column  which  sometimes  may  exceed  safe  levels  as  is  evidenced  by  the 
large  number  of  vertebral  fractures  which  have  been  reported  following 
ejection. 

Moffatt  and  Howard  (1968)  made  an  extensive  compilation  of  U.S. 
Air  Force  and  Navy  aircraft  ejections  during  the  period  1959-1967. 

They  found  that  17%  of  the  ejections  resulted  in  vertebral  fractures 
with  most  (70%)  occurring  in  the  lower  thoracic  region  (T7  through 
T12)  (Figure  2).  Shannon  (1970)  studied  561  combatant  and  noncom¬ 
batant  ejections  which  occurred  during  1967  and  1968.  In  this  study, 
8.5%  of  the  ejectees  demonstrated  vertebral  fractures  which  were 
caused  entirely  by  the  ejection  system.  Again  the  lower  thoracic  and 
upper  lumbar  regions  were  found  to  be  the  prime  injury  sites. 

At  a  recent  (1973)  Ejection  Seat  Committee  meeting  of  AGARD,  the 
continuing  seriousness  of  the  ejection  problem  was  reported  by  Prof. 

R.  P.  Delahaye  of  France  on  behalf  of  all  the  participating  NATO  coun¬ 
tries.  These  included  the  British,  Hellenic,  Italian,  German  and  French 
Air  Forces.  The  U.S.  was  represented  by  statistics  from  the  Army  and 
Air  Force,  but  not  the  Navy.  Of  678  ejection  episodes  in  1972,  there 
were  114  deaths;  98  pilots  sustained  spinal  injuries,  which  resulted  in 
160  vertebral  fractures.  The  distribution  of  160  fractures  in  98  pilots 
is  shown  in  Figure  3.  Figure  4  shows  the  distribution  of  fractures  in 
pilots  with  multiple  spinal  fractures.  Figures  3  and  4  show  the  great¬ 
est  incidence  of  injury  to  be  at  T7  -  T8  and  at  T12  -  LI,  the  thoraco¬ 
lumbar  transition  vertebra.  Prof.  Delahaye  further  stated  that  the  fre¬ 
quency  of  multiple  fractures  was  increasing.  There  exists  a  general 
consensus  that  poor  initial  configuration  is  the  principle  cause  of 
fracture,  and  that  this  has  been  proven  convincingly  through  the  use  of 
the  powered  inertial  reel  in  the  German  Air  Force.  With  respect  to  an¬ 
thropomorphic  type,  he  reported  that  "thin  and  tall"  pilots  have  higher 
frequency  of  fracture  in  the  T7  -  T8  region,  while  "stout  and  medium" 
pilots  tend  to  have  T12  -  LI  fractures. 

In  addition  to  the  medical  aspects,  there  is  a  substantial  eco¬ 
nomic  cost  associated  with  unnecessary  fracture.  Ewing  (1971)  calcu¬ 
lated  an  average  yearly  cost  to  the  U.S.  Navy  of  $6,797,718  for  aviators 
who  had  sustained  vertebral  fractures. 

Latham  (1957)  was  one  of  the  first  investigators  to  correlate 
vertebral  fracture  with  the  biomechanics  of  the  spinal  column.  He 
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Figure  1.  Martin-Baker  Acceleration  Pulse 
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Figure  2.  Spinal  Injury  Statistics: 

Moffatt  and  Howard 
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Figure  3.  Spinal  Injury  Statistics:  Delahaye 
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reported  on  work  done  by  the  German  Air  Force  during  the  early  40' s 
which  ascribed  anterior  fractures  of  the  vertebrae  to  excessive  flexion 
of  the  trunk.  Additionaly,  he  noted  the  need  for  a  lumbar  pad  to  at¬ 
tain  correct  spinal  alignment  prior  to  ejection  and  the  increased  bend¬ 
ing  induced  by  the  rearward  rotation  of  the  thrust  axis  relative  to  the 
geometric  axis  of  the  spine.  Bosee  and  Payne  (1961)  studied  the  Martin- 
Baker  ejection  system.  They  postulated  that  the  common  anterior  lip 
fracture  was  due  to  bending  stresses  which  could  be  induced  by  poor  in¬ 
itial  spinal  alignment.  This  hypothesis  was  later  corroborated  experi¬ 
mentally  by  Vulcan  et  al.  (1970).  They  attached  strain  gages  to  cada¬ 
ver  vertebrae  and  demonstrated  the  existence  of  large  compressive 
stresses  on  the  anterior  edge  of  the  vertebral  body. 

Mathematical  models  have  also  been  employed  to  study  the  verte¬ 
bral  fracture  problem.  Early  models  were  derived  to  study  only  axial 
deformation.  Latham  (1957)  modelled  the  man,  ejection  seat,  spine  con¬ 
figuration  as  a  single  spring-mass  system.  Unfortunately,  the  use  of  a 
single  spring  for  the  spine  prevents  any  study  of  regional  spinal 
failures,  e.g.,  the  propensity  of  lower  thoracic  and  upper  lumbar  frac¬ 
ture  noted  earlier.  This  deficiency  was  somewhat  rectified  by  Toth 
(1967)  who  used  eight  spring,  mass  elements  to  study  the  Til-pelvis  seg¬ 
ment.  However,  the  model  was  still  confined  to  a  study  of  axial  de¬ 
formations. 

Vulcan  and  King  (1970)  used  a  three  mass  model  to  study  the  ef¬ 
fect  of  harness  and  seat  back  constraints  (Figure  5).  The  force  levels 
at  the  lower  vertebral  levels  were  determined  by  changing  the  dimen¬ 
sions  and  mass  of  the  torso  unit  to  correspond  to  different  spinal 
levels.  Again  it  was  found  that  large  bending  and  axial  forces  existed 
in  the  lower  thoracic  and  upper  lumbar  spine. 

Ome  and  Liu  (1971)  have  derived  the  most  comprehensive  model  of 
the  spinal  column  presently  available.  This  model  consists  of  a  series 
of  rigid  masses  connected  by  viscoelastic,  massless  beams.  This  formu¬ 
lation  rectified  most  of  the  earlier  shortcomings  by  accounting  for  axi¬ 
al,  shear,  and  bending  deformations  of  the  spinal  column,  the  variable 
size  of  vertebrae  and  discs,  the  natural  curvature  of  the  spine,  and  the 
eccentric  inertial  loading  on  the  spine  produced  by  the  mass  of  the  head 
and  trunk.  A  study  of  a  typical  unrestrained  ejection  sequence  with 
this  model  indicated  that  large  anterior  stresses  did  exist  in  the  lower 
thoracic  and  upper  lumbar  regions  if  the  stresses  due  to  axial  and  bend¬ 
ing  forces  were  combined  according  to  the  strength  of  materials  formula, 
cr  =  P/A  +  Mc/I.  A  serious  shortcoming  of  this  model,  however,  was  the 
fact  that  it  did  not  account  for  the  restraining  effect  of  the  shoulder 
harness  and/or  seat  back. 

The  purpose  of  the  present  investigation  was  to  extend  the  Ome- 
Liu  mod, el  by  including  the  effects  of  external  constraints.  This  was 
accomplished  and  the  model  used  to  study  the  effects  of  initial  curva¬ 
ture  on  spinal  injury  potential  of  the  MK-J5(D)  ejection  system. 
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Figure  5.  Simplified  Spinal  Model 

(Taken  from  Vulcan  §  King) 
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MK-J5  EJECTION  SYSTEM 


The  separation  of  the  aviator  from  the  cockpit  is  accomplished 
by  the  Model  MK-J5  ejection  seat  shown  in  Figure  6.  The  MK-J5  is  a 
fully  automatic,  cartridge  operated  ejection  system  which  incorporates 
the  following  features: 

1.  A  high  velocity  ejection  gun  which  provides  sufficient  accelera¬ 
tion  to  clear  the  vertical  stabilizer  and  attain  sufficient  al¬ 
titude  for  safe  chute  deployment  during  low  level  ejection. 

2.  A  dual  automatic  leg  restraining  mechanism  which  helps  prevent 
flail  injuries  to  the  leg. 

3.  An  integrated  back  pad  and  seat  pan. 

4.  A  face  blind  which  helps  restrain  forward  head  motion. 

Despite  these  features,  several  shortcomings  of  the  system  have 
been  described.  Bosee  and  Payne  (1961)  and  Kaplan  (1972)  have  indicated 
that  the  MK-J5(A,B)  ejection  system  forced  the  aviator  into  a  flexed  po¬ 
sition  in  the  lower  spinal  region  and  also  placed  the  mass  of  the  head 
and  thorax  further  anterior  than  would  be  commensurate  with  good  spinal 
alignment.  This  initial  configuration  would  predispose  the  aviator  to 
spinal  fracture.  Additionally,  Kaplan  (1972)  indicated  that  insuffici¬ 
ent  thigh  support  could  exacerbate  flexion  in  the  lumbar  region  through 
an  induced  rotation  of  the  pelvis.  He  also  pointed  out  the  necessity 
for  proper  design  of  the  lumbar  support  pad. 

In  order  to  rectify  some  of  these  problems,  modifications  were 
made  to  the  MK-J5(A,B)  system.  Kaplan  (1972)  reported  on  these  changes 
which  are  designated  the  MK-J5(D)  system  (Figure  7).  It  was  found  that 
comfort  had  increased  through  the  incorporation  of  an  improved  harness 
system  and  an  elongation  of  the  thigh  support.  In  addition,  thoracic 
flexion  was  reduced  when  compared  to  that  of  the  earlier  MK-J5(A,B)  sys¬ 
tem,  but  certain  size  percentile  aviators  were  still  predisposed  to 
fracture  due  to  improper  alignment.  Kaplan  (1972)  notes:  "The  5th  and 
95th  percentiles  (sitting  height)  appear  to  be  more  predisposed  to  ver¬ 
tebral  fracture  than  the  40th  percentile  due  to  seat  back  contour  de¬ 
sign.  The  face  curtain  firing  position  increased  vertebral  flexion  in 
the  95th  percentile." 

Kaplan  (1972)  also  analyzed  statically  the  probability  of  frac¬ 
ture  from  spinal  alignment  x-rays  taken  while  5th,  40th  and  95th  per¬ 
centile  (sitting  height)  aviators  were  seated  in  the  MK-J5 (D)  system  in 
the  face  curtain  (position  1)  or  secondary  "D"  ring  position  (position 
2).  A  consideration  of  the  spinal  alignments  indicated  that  the 
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Figure  6.  MK-J5(A,B)  Ejection  System 
(Taken  from  Kaplan) 
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Figure  7.  MK-J5(D)  Ejection  System 
(Taken  from  Kaplan) 
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probable  level  of  fracture  differed  with  percentile  sitting  height  as 
given  below  in  Table  1: 


TABLE  1 


PREDICTED  FRACTURE 

LEVELS,  Kaplan 

(1972) 

Percentile  Sitting 

Height 

5^ 

40 

95 

Spinal  Level  T12 

T9- 11 

T9 

In  order  to  study  these  statically  derived  observations'  under  dynamic 
conditions,  the  Ome-Liu  model  of  the  spinal  column,  with  restraint 
added,  was  used  to  investigate  the  dynamic  response  of  the  spine  using 
data  reported  by  Kaplan  (1972). 

ORNE-LIU  MODEb  WITH  CONSTRAINTS 

The  main  load  carrying  component  of  the  human  body  during  ejec¬ 
tion  is  the  spinal  column  (Figure  8) .  It  consists  of  33  vertebrae  and 
23  intervertebral  discs.  In  accordance  with  their  weight  bearing  func¬ 
tions,  the  vertebrae  and  discs  increase  in  size  caudally.  Regional 
variations  in  the  shape  of  these  components  account  for  the  cervical, 
thoracic  and  lumbar  curvatures  found  in  the  adult  spine.  The  mobility 
of  the  spine  is  primarily  due  to  the  flexibility  of  the  intervertebral 
discs  while  the  limits  of  motion  are  defined  by  ligaments,  facets,  mus¬ 
cles  and  bony  overlaps.  An  important  consideration  in  the  present  in¬ 
vestigation  is  the  effect  of  the  ribs  in  preventing  midsagittal  flexion 
of  the  thoracic  column.  The  relatively  stiff,  box-like  thoracic  cage 
has  practically  a  rigid-body  motion  relative  to  the  cervical  or  lumbar 
spine,  i.e.,  the  mobility  of  the  spine  in  the  midsagittal  plane  is 
largely  confined  to  the  lumbar  and  cervical  regions. 

Basing  their  work  on  the  anatomy  of  the  spine,  Orne  and  Liu  (1971) 
have  formulated  the  most  advanced  dynamic  spinal  model  now  available. 

The  model  consists  of  alternating  rigid  masses,  which  represent  the  in¬ 
ertial  properties  of  the  torso,  and  massless  deformable  links  which  ac¬ 
count  for  the  elastic  properties  of  the  spinal  column  (Figure  9).  This 
discretization  allows  the  spinal  model  to  account  for,  simultaneously: 

1.  the  axial,  shear  and  bending  deformations  of  the  discs, 

2.  the  variable  size  of  vertebrae  and  discs, 

3.  the  natural  curvature  of  the  spine, 

4.  the  eccentric  inertial  loading  on  the  spine  produced  by  the 
mass  of  the  head  and  trunk. 
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.  Figure  9.  Ome-Liu  Discrete  Parameter  Model  of  the  Spinal  Column 


Figure  10.  Free-Body  Diagram  of  Vertebral  Body 


'> 
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The  following  equations  of  motion  for  the  model  are  given  only 
for  the  purpose  of  showing  how  the  constraints  are  added.  Details  of 
the  derivation  may  be  found  in  Ome  and  Liu  (1971).  The  forces  acting 
on  the  vertebral  body  are  shown  in  Figure  10,  which  is  the  free-body 
diagram  for  a  typical  vertebral  segment.  The  center  of  mass,  mi,  has 
an  eccentricity,  ei,  with  respect  to  the  vertical  axis  of  the  vertebral 
body.  Ji  is  the  mass  moment  of  inertia  about  an  axis  perpendicular  to 
the  sagittal  plane  and  through  the  center  of  mass  and  hi  is  one-half 
the  vertebral  height.  The  equations  of  motion  of  the  vertebrae  are: 


and 


and 


m.  x.  =  Q.  -  Q. 

1  1  xi  \L-1 


m.  z.  =  P.  -  P.  , 
li  i  l-l 


where 


J. 

1 

<t>.  =  (F3.  ,  • 
l  v  l-l 

-  Bi> 

-  C-.Q.  .  +  C-.Q.  +  C, .P.  .  -  C..P. 

2ixi-l  3ixi  li  l-l  4i  l 

cli 

=  -  (h.  sin 

l 

*i  + 

e.  cos  4.) 
l  Ti 

C2i 

=  -  (h.  cos 

l 

♦i  - 

e.  sin  ) 

l  yi 

C  „ . 

=  h.  cos  d>. 

+  e. 

sin  4. 

3i 

i  Ti 

i 

Ti 

C.. 

=  h.  sin  4. 

-  e. 

cos  4. 

4i 

i  l 

i 

Ti 

The  Qi 

,  Pi,  F3i,  and  Bi 

are  the  force  and  moment  reactions  of  the 

discs  on  the  vertebrae.  To  these  equations  must  be  added  the  constrain¬ 
ing  forces  of  the  shoulder  harness,  face  curtain,  and  seat  back.  Al¬ 
though  the  Orne-Liu  model  can  be  used  to  represent  each  vertebral  level 
as  a  separate  entity,  it  would  unnecessarily  complicate  the  modelling 
process  to  discretize  the  entire  column.  Thus,  the  segments  T1  to  T6 
were  modelled  as  a  rigid  link.  This  assumption  is  based  on  the  work  of 
Snyder  et  al.  (1971),  who  found  that  the  thoracic  segments  which  support 
the  rib  cage  act  as  a  rigid  link  insofar  as  midsagittal  bending  is  con¬ 
cerned.  Also,  the  vast  majority  of  vertebral  fractures  occur  below  the 
level  of  T6.  Additionally,  the  cervical  region  was  lumped  into  two  seg¬ 
ments.  One  segment  contained  vertebrae  C7  through  C4  and  the  other  C3 
through  Cl .  Linear  springs  were  used  to  model  the  restraining  effects 
of  the  shoulder  harness,  face  curtain,  and  seat  back.  By  including  the 
restraining  forces  the  equations  of  motion  change  to: 


mx.  =  Q.  -  Q.  -  RX. 

i  xi  Hi-1  l 


mz.  =  P.  -  P.  ,  -  RZ. 
l  i  l-l  l 
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where 


J.  X.  =  (F3.  .  -  B.  -  MX.  -  MZ.)  -  CL.  Q.  ,  +  CL.  Q. 

1  Yi  *■  l-l  1  l  xJ  2i  xi-l  3i 


+  C. .  P.  .  -  CL.  P. 
li  l-l  4i  i 


RX.  =  k  (u.  -  u.  +  d.  sin  A.  -  d.  sin  A .  } 
1  xv  1  10  1  Y1  1  Yio 

RZ.  =  k  (w.  -  w.  +  d.  cos  (j>.  -  d.  cos  <j>.  ) 
1  zv  1  lO  1  1  1  10 


HXi  =  -  RXi  (Cn  *  di  cos  <fr.) 


MZ.  =  RZ.  (CL.  +  d.  sin  <f>.) 
l  l  *■  li  l 


The  quantities,  RXi,  RZi,  MXi  and  MZi  are  the  component  force 
and  moment  reactions  of  the  restraining  springs  on  the  vertebrae,  di  is 
the  distance  along  the  vertebral  axis  from  the  inferior  vertebral  end 
plate  to  the  point  of  attachment  of  the  restraining  spring  and  ui  and 
wi  are  the  instantaneous  configuration  co-ordinates.  For  the  seat  back 
and  face  curtain  restraints,  di  was  taken  as  hi,  one-half  the  vertebral 
height.  The  shoulder  harness  springs  were  attached  at  the  upper  and 
lower  ends  of  the  T1  through  T6  segment.  Thus  di  was  set  to  2hi  for 
the  upper  harness  spring  and  zero  for  the  lower.  The  spring  constants 
kx  and  kz  were  selected  on  the  basis  of  experimental  work  done  by  Vulcan 
and  King  (1970).  Their  tests  had  indicated  an  effective  kx  of  200  lb/in 
when  the  torso  displacement  and  horizontal  shoulder  strap  force  were 
considered.  The  vertical  spring  constant  kz  was  chosen  as  50  lb/in 
since  Vulcan  and  King  (1970)  noted  an  approximately  150  lb  decrease  in 
the  vertical  portion  of  the  shoulder  strap.  In  the  absence  of  informa¬ 
tion  on  the  effective  elasticity  of  the  seat  back  and  face  curtain  re¬ 
straint,  the  shoulder  strap  spring  constants  kz  and  kx  were  utilized 
for  these  restraining  springs  as  well.  These  constants  appeared  to 
give  reasonable  results  when  used  in  the  face  curtain  firing  position. 
However,  when  the  5th  percentile,  secondary  "D"  ring  configuration  was 
analyzed,  excessive  forward  displacement  and  rotation  of  the  torso  took 
place  due  to  the  large  forces  generated  by  the  free  head.  Several 
spring  constants  were  utilized  for  the  shoulder  strap  until  one  which 
would  allow  only  a  relatively  small  lateral  displacement  of  the  torso 
(one  or  two  inches)  was  found.  The  spring  constant,  kx,  utilized  for 
the  shoulder  strap  with  the  head  free  was  1400  lb/in. 


One  additional  restraint  was  required  to  account  for  the  forces 
developed  when  the  head  strikes  the  anterior  chest  wall.  This  restraint 
was  modelled  as  a  linear  spring  which  was  engaged  when  the  head  rotated 
45  degrees  forward  relative  to  the  torso.  The  spring  constant  was  se¬ 
lected  so  that  the  head  would  not  rotate  more  than  10  degrees  further 
after  engaging  the  spring.  The  spring  constant  so  chosen  was  100  lb/ra- 
dian. 

With  these  modified  equations  of  motion,  the  Ome-Liu  model  was 
used  to  study  the  MK-J5(D)  ejection  system. 
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BIOMECHANICAL  DATA 


Prior  to  exercising  the  constrained  Ome-Liu  model  of  the  spine, 
the  physical  parameters  of  the  model  should  be  specified.  For  the  pur¬ 
pose  of  this  investigation,  the  material  properties  assumed  by  Ome  and 
Liu  (1971)  for  shear,  axial  and  bending  behavior  of  the  disc  will  be 
used.  Thus,  the  viscoelastic  behavior  of  the  intervertebral  disc  is 
governed  by  a  single  set  of  material  parameters:  p^,  q0,  q^,  which 
have  the  values  0.030  sec,  1100  psi,  and  72  psi-sec  respectively.  The 
values  of  G,  the  shear  modulus,  and  E^,  the  bending  modulus  were  taken 
as  2200  psi  and  7700  psi  respectively.  The  shape  factor,  F,  normally 
associated  with  shear,  was  taken  to  be  1.333  (circular  cross-section 
assumed) . 

The  remaining  input  data  was  taken  and/or  derived  from  the  re¬ 
sults  of  Kaplan  (1972).  The  initial  values  of  spinal  curvature,  <^0, 
were  taken  from  the  spinal  alignments  of  the  aviators  as  seated  in  the 
MK-J5(D)  seat  (Figure  11).  It  should  be  noted  that  the  angles  reported 
by  Kaplan  (1972)  were  taken  with  respect  to  the  thrust  axis  and  not  the 
vertical  axis.  For  the  MK-J5  seat,  the  thrust  axis  is  rotated  rearward 
19  degrees  from  the  vertical.  Also,  for  one  case  (5th  percentile)  the 
lumbar  alignments  for  an  improved  lumbar  support  pad  described  by 
Sanford  and  Kellet  (1966)  were  utilized.  The  vertebral  height  was 
taken  as  the  distance  between  the  centers  of  the  vertebral  bodies  as 
measured  by  Kaplan  (1972),  minus  the  disc  height  used  by  Ome  and  Liu 
(1971).  The  calculation  of  mi,  ei,  and  Ji  was  somewhat  more  complicated. 
Liu  and  Wickstrom  (1972)  reported  regression  equations  for  mi,  ei,  and 
Ji  derived  from  segmented  cadavers.  These  equations  give  the  inertial 
property  distribution  of  the  human  torso  in  a  form  that  can  be  used  in 
the  Ome-Liu  model.  However,  in  order  to  use  the  regression  equations, 
the  height,  weight,  breadth  and  chest  depth  of  the  individual  are  re¬ 
quired.  Kaplan  (1972)  reported  measurements  of  height,  weight  and  sit¬ 
ting  height.  The  additional  required  measures  were  determined  using  the 
percentile  data  given  by  Kaplan  (1972)  and  anthropometric  measurements 
collected  by  Schane  et  al.  (1969).  The  final  regression  parameters  are 
given  in  Table  2. 

TABLE  2 


INERTIAL  DISTRIBUTION  REGRESSION  PARAMETERS 


Percentile 
[Sitting  Height] 

Weight 

[lbs] 

Height 

[in] 

Chest  Depth 
[in] 

Hip  Breadth 
[in] 

5 

144.1 

64.3 

20.0 

31.83 

.  40 

172.0 

67.9 

21.1 

33.23 

95 

189.5 

75.5 

22.1 

34.63 

While  these  parameters  define  mi,  e^,  and  Ji  for  the  torso,  it  is  not 
possible  to  calculate  the  cervical  or  head  inertial  properties  from  the 
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MCE  CURTAIN  SECONDARY  “D"EINO 


Figure  11.  Spinal  Alignment,  40th  Percentile 
Aviator  (Taken  from  Kaplan) 
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regression  equations.  Thus,  the  cervical  and  head  inertial  data  are 
taken  from  Pontius  et  al.  (1972).  An  additional  consideration  in  the 
present  investigation  is  the  mass  of  the  helmet  which  is  worn  by  the 
aviator.  The  physical  properties  associated  with  the  helmet  were  de¬ 
termined  by  the  methods  of  Walker  et  al.  (1973)  and  Liu  et  al.  (1971). 

To  find  the  center  of  gravity  (c.g.)  of  the  helmet,  with  visor  down,  it 
was  hung  from  two  points  in  the  midsagittal  plane  (Figure  12) .  It  was 
seen  that  the  c.g.  of  the  helmet  approximated  the  location  of  the  c.g. 
of  the  head  as  given  by  Walker  et  al.  (1973).  Thus,  it  was  assumed 
that  the  helmet  mass  (0.088  lb-sec2/in)  and  mass  moment  of  inertia 
(0.118  lb-sec2/in)  could  simply  be  added  to  the  respective  values  for 
the  head. 

Using  the  regression  parameter  values  given  in  Table  2,  the  re¬ 
gression  equations  of  Liu  and  Wickstrom  (1972),  and  the  data  given  in 
Kaplan  (1972),  the  complete  set  of  input  data  for  the  5th,  40th,  and 
95th  percentile  aviators  were  determined  (Tables  3-5). 

RESULTS  AND  DISCUSSION 

The  idealized  trapezoidal  acceleration  pulse  applied  vertically 
at  the  base  of  the  spine  is  shown  in  Figure  13  along  with  an  actual  re¬ 
cord  of  a  MK-J5(D)  acceleration  time  history.  The  trapezoidal  accelera¬ 
tion  pulse  was  used  to  study  several  configurations  of  the  model  as 
given  in  Table  6.  The  configurations  were  selected  mainly  on  the  basis 
of  spinal  alignment  data  reported  by  Kaplan  (1972)  for  the  face  curtain 
(position  1)  and  secondary  "D"  ring  (position  2)  ejection  positions  of 
three  aviators  seated  in  the  MK-J5(D)  ejection  seat.  Three  aviators, 
5th,  40th  and  95th  percentile  men  in  terms  of  their  sitting  height,  were 
selected  to  represent  their  anthropometric  class.  The  constraints  were 
varied  so  that  the  effects  of  seat  back,  face  curtain  and  shoulder  har¬ 
ness  restraints  could  be  studied.  Additionally,  the  lumbar  spinal 
alignment  induced  by  an  improved  McDonnell  lumbar  pad  discussed  by 
Sanford  and  Kellet  (1966)  was  studied. 

The  first  configuration  served  to  define  the  basic  unconstrained 
response  of  the  spinal  column  to  caudocephalad  acceleration  and  was 
used  as  a  baseline  to  observe  changes  induced  by  the  addition  of  re¬ 
straints.  The  graphical  results  shown  in  Figures  14  through  18  give 
the  unconstrained  response  of  the  5th  percentile  aviator.  In  these  and 
all  following  figures,  axial  compressive  forces,  bending  moments  which 
induce  compression  in  the  anterior  half  of  the  vertebral  body,  and 
shearing  forces  directed  anteriorly  are  plotted  with  positive  ordinates. 
As  could  be  expected  of  the  case  without  restraint,  the  torso  flexes 
rapidly  and  reaches  a  horizontal  position  at  170  ms.  In  order  to  study 
the  regional  variation  in  forces,  an  "effective"  normal  compressive 
stress  was  calculated  according  to  the  strength  of  materials  formula: 
ac  =  P/A  ±  Mc/I.  Positive  values  of  this  function  indicate  the  exis¬ 
tence  of  net  compressive  stress  on  the  anterior  half  of  the  vertebral 
body.  Additionally,  an  effective  axial  stress,  a a  =  P/A,  was  calcula¬ 
ted.  Both  effective  stress  functions  were  calculated  using  the  forces 
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(b) 

Figure  12.  Center  of  Gravity  of  Aviator's  Helmet 
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TABLE  3 


5TH  PERCENTILE  ANTHROPOMETRIC  DAT  A 


*-******•******************* ************************************************************************* •••*•••*•*** 

MCDONNELL 


VERTEBRAL 

level 

DISC 
AREA 
(  IN1) 

f* ******** 

DISC 
HEIGHT 
(  IN  ) 

«*  ******** 

area  moment 
of  inertia 
(  IN1* ) 

************< 

MASS 

(LB- SEC*  /IN) 

i***********i*4 

ROTATORY 

inertia 

(LB-SEC^-IN) 

.************ 

vertebral 

HEIGHT 

(IN) 

**********4 

ECCENTRICITY 

(IN) 

!>****♦♦***♦*** 

POS  1 
ANGIE 
(OEG) 

80S  2 
angle 
( DEI  1 

LUMBAR  PAD 

angle 

( OEG  1 

L6 

2.72 

0 .36 

0.306 

0.0118 

0. 0528 

1.32 

0.89 

-4. 

-6. 

23. 

L4 

2.66 

0.48 

0.220 

0.0106 

0 . 04B 1 

0.98 

1  .06 

-8. 

-14. 

24. 

L3 

2.44 

0.42 

0.282 

0.0104 

0.0472 

1.29 

1.19 

-5. 

-6. 

22. 

L2 

2.22 

0.41 

0.  169 

0.0091 

0.0442 

1.12 

1.11 

-2. 

-6. 

21. 

Li 

1.94 

0.38 

0. 135 

O.0092 

0.0518 

1.02 

1.17 

2. 

0. 

19*- 

T12 

1  .97 

0.28 

0.119 

0.0083 

0. 0490 

1.03 

1.37 

7. 

6. 

1  7. 

Til 

1.84 

0.17 

0.  U5 

0.0067 

0.0383 

1.13 

1.36 

7. 

9. 

15. 

no 

1 .52 

0.17 

0. 102 

0.0057 

0.0326 

1.15 

1  .44 

18. 

16. 

13. 

T9 

1.47 

0.15 

0.097 

0.0063 

0.0365 

0.85 

1.51 

17. 

U. 

11. 

T  8 

1  .  36 

0.15 

C  .  DSC 

0.0956 

0.0339 

0.85 

1.26 

8. 

19. 

8. 

T? 

1.32 

0. 15 

0.065 

0.0050 

0.0320 

0.76 

1  .08 

23. 

19. 

16. 

T61 

1 . 30 

0.14 

0.  060 

0.03?4 

0.2246 

5.  16 

1.10 

23. 

45. 

23. 

C  74 

1.20 

0.24 

0.046 

0.0152 

0.0564 

3.11 

0.00 

23. 

10. 

23. 

C  3 1 

1.12 

0.24 

0.036 

0.0030 

0.0064 

1.26 

0.00 

20. 

7  . 

20. 

H 

0.49 

0.2  4 

0.013 

0.0426 

0.51 30 

1  .  10 

0.70 

0. 

0. 

0. 

40TH  PERCFNTILE  ANTHROPOMETRIC  DATA 


**********1******************************  **********  ***********!*f*************************5!!*:****':******** 
VERTEBRAL  OISC  DISC  AREA  MOMENT  ROTATORY  VERTEBRAL  POS  I 

LEVEL  AREA  HEIGHT  OF  INERTIA  MASS  INERTIA  HEIGHT  ECCENTRICITY  ANGLE 

(IN4)  (IN)  (IN*)  (LB-SEC2/IN)  (LB-SEC-IN)  (IN)  (IN)  (OEG) 


************************************************ ***************** ****** ****************************** 


L5 

2.72 

0.36 

0.306 

0.0132 

0.0720 

1.32 

0.33 

-1  . 

L4 

2.66 

0.48 

0.220 

0.0125 

0.0698 

1.38 

0.59 

-6, 

L3 

2.44 

0.42 

0.282 

0.0121 

0.0689 

1.29 

1.06 

0. 

L2 

2.22 

0.41 

0.169 

0.0110 

0.0641 

1.12 

0.93 

-3. 

LL 

1.94 

0.38 

0.135 

0.0107 

0.0669 

1.22 

1  .04 

-3. 

T12 

1.97 

0.28 

0.119 

0.0093 

0.0599 

1.33 

1.20 

-11  . 

Til 

1.84 

0.17 

0.115 

0.0086 

0.0574 

0.83 

1.21 

-10. 

TIC 

1.52 

0.17 

0.  102 

0.0071 

0.0484 

1.15 

1.38 

0. 

T9 

1.47 

0.15 

0.097 

0.0077 

0.0522 

1.10 

1.44 

10  ■, 

T8 

1.36 

0.15 

0.080 

0.0072 

0.0489 

0.95 

1.13 

9. 

T7 

1.32 

0.15 

0.065 

0.0071 

0.0515 

0.76 

1.01 

6. 

T6i 

1.30 

0.14 

0.060 

0.0359 

0.2558 

4.51 

1.05 

27. 

C74 

1.20 

0.24 

0.046 

0.0152 

0.0553 

2.36 

0.00 

21  . 

C31 

1.12 

0.24 

0.038 

0.0030 

0.0Q63 

0.96 

o.oc 

11. 

H 

C  .  49 

0.24 

0.013 

0.0426 

0.5130 

1.10 

0.70 

0. 

TABLE  5 


95TH  3ERCENTILE  ANTHROPOMETRIC  DATA 


Jfcjjr  *  +  *:$$:*:£*  ^JiifcfcJt;**^*:**:*#*#***#**  £  *  $  $  **  *  <=  fc  *  *  >r  <=  $$#$*******  ’ft*#**#***  A*  *£«"**=*  *$**£*****  $$#*=!■:{: 

VERTEBRAL  DISC  DISC  AREA  MOMENT  ROTATORY  VERTEBRAL  POS  1 

LEVEL  AREA  HEIGHT  OF  INERTIA  MASS  INERTIA  HEIGHT  ECCENTRICITY  ANGLE 

(  IN'2- )  (IN)  (IN*)  <LB-SEC2/IN)  (LB-SEC^-IN)  (IN)  (IN)  (DEG) 

fc*******:#*#*  *#*««*<:*<!«**  $  $  *  *  *  $  fr  #  *  *  *  *  $  #  ♦♦  ♦  $  *  #  *  fc  $  *  #  *  ^  $  **  ♦  *  ♦  £ 


L5 

2.72 

0.36 

0.306 

0.0152 

0.0879 

1.32 

0.72 

11  . 

L4 

2.66 

0.48 

0.220 

0.0153 

0.0878 

1.38 

0.92 

-1  . 

L3' 

2.44 

0.42 

0.282 

0.0146 

0.06  70' 

1.39 

1.15 

-3. 

L2 

2.22 

0.41 

0.169 

0.0138 

0.0807 

1.32 

1.07 

-13. 

Ll 

1.94 

0.38 

0.  135 

0.0130 

0.0795 

1.42 

1.13 

-3. 

T12 

1.97 

•  0.28 

0.119 

0.0109 

0.0691 

1.43 

1.32 

-11. 

TIL 

1.84 

0.17 

0.115 

0.0115 

0.0733 

1.08 

1.31 

-18. 

T10 

1.52 

0.17 

0. 102 

0.0092 

0.0615 

1.15 

1.42 

-15. 

T9 

1.47 

0.15 

0.09  7 

0.0098 

0. 0652 

0.95 

1.49 

-10. 

TB 

1.36 

0. 15 

0.080 

0.0096 

0.0615 

0.95 

1.22 

0. 

T7 

1.32 

0.15 

0.065 

0.0103 

0.0678 

0.96 

1.05 

5. 

T61 

1.30 

0.14 

0. 060 

0.0412 

0.3519 

5.26 

1.09 

5. 

C74 

1.20 

0.24 

0.046 

0.0152 

0.0577 

3.01 

0.00 

19. 

C3L 

1.12 

0.24 

0.038 

0.0030 

0.0064 

1.  16 

0.00 

-2  . 

H 

0.49 

0.24 

o . :  1 3 

0.0426 
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Figure  13.  MK-J5(D)  Acceleration  Pulse 
with  Trapezoidal  Idealization 
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Figure  14.  Spinal  Alignment,  Sth  Percentile, 
Position  1,  Unrestrained 
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TABLE  6 


MODEL  CONFIGURATIONS 


CONFIGURATION 

PERCENTILE  MAN 
(Sitting  Height) 

INITIAL 

SPINAL 

POSITION 

RESTRAINTS 

1 

5 

1 

None 

2(a)i 

5 

1 

Face  Curtain,  Shoulder 
Harness 

2(b) 

40 

1 

Face  Curtain,  Shoulder 
Harness 

2(c) 

95 

1 

Face  Curtain,  Shoulder 
Harness 

3(a) 

5 

1 

Face  Curtain,  Shoulder 
Harness,  Seat  Back 

3(b) 

40 

1 

Face  Curtain,  Shoulder 
Harness,  Seat  Back 

3(c) 

95 

1 

Face  Curtain,  Shoulder 
Harness,  Seat  Back 

4 

5 

2 

Shoulder  Harness,  Seat  Back 

5 

5 

McDonnell 

Face  Curtain,  Shoulder 

Pad 

Harness,  Seat  Back 
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and  moments  existing  on  the  superior  surface  of  the  rigid  vertebral 
body.  Figure  IS  shows  ac  for  the  unconstrained  5th  percentile  aviator. 
The  maximum  normal  stress  is  seen  to  occur  at  L5  but  this  peak  is  due  to 
the  clamped  beam  boundary  condition.  Thus  the  secondary  peak  at  L2  is 
more  indicative  of  the  location  of  maximum  stress.  Figure  16  shows  the 
variation  of  maximum  axial  force  and  effective  axial  stress  for  the  un¬ 
constrained  case.  The  maximum  axial  force  is  found  at  the  base  of  the 
spine  but  the  maximum  axial  stress  occurs  in  the  lower  and  middle  thor¬ 
acic  region.  This  was  found  to  be  true  for  all  cases  studied.  In 
several  instances  (Table  7) ,  maximum  axial  stresses  did  occur  in  the 
upper  lumbar  region  but  when  this  was  the  case  secondary  maximums,  usu¬ 
ally  only  several  percent  less  in  magnitude,  wepe  found  higher  in  the 
column.  The  axial  stress  was  relatively  constant  along  the  column.  It 
appears  that  the  variation  of  cross-sectional  area  of  the  discs  and  ver¬ 
tebrae  is  such  that  a  constant  axial  stress  is  maintained  in  the  column. 
The  time  histories  of  the  axial  force,  bending  moment,  and  shear  force 
on  the  superior  surface  of  LI  are  shown  in  Figure  17.  The  time  of  maxi¬ 
mum  normal  stress  corresponds  to  the  occurence,  of  peak  bending  moment 
since  the  axial  force  is  relatively  constant.  Figure  18  shows  the  axial 
force  and  bending  moment  distribution  along  the  spinal  column  at  170  ms. 

Results  for  Configuration  2  are  given  in  Figures  19  through  23. 
Figure  19  shows  the  displacements  of  the  spinal  column  when  the  lower 
vertebrae  are  not  constrained.  Comparing  Figures  14  and  19,  we  see  that 
the  shoulder  harness  and  face  curtain  prevent  flexion  of  the  torso. 

Thus,  the  resulting  deformation  of  the  column  is  restricted  mainly  to 
the  lower  thoracic  and  lumbar  regions.  This  serves  to  place  the  maximum 
normal  stress  at  points  of  maximum  spinal  curvature  as  can  be  seen  from 
Figure  20.  It  should  be  emphasized  at  this  point  that  the  displacement 
curve  is  the  loci  of  all  the  midpoints  of  the  inferior  vertebral  end- 
plates  while  the  forces  are  taken  with  respect  to  the  superior  endplates. 
The  5th  percentile  aviator  has  a  maximum  normal  stress  at  LI,  the  40th 
at  T12 ,  and  the  95th  at  T9.  These  dynamic  stress  levels  coincide  al¬ 
most  exactly  with  the  static  observations  of  Kaplan  (1972)  who  predicted 
that  maximum  stress  levels  would  occur  at  T12,  Til,  and  T9,  respectively. 
Thus,  the  curvatures  as  determined  from  spinal  alignment  x-rays  under 
static  conditions  are  useful  for  determining  spinal  injury  levels  for 
this  ejection  configuration.  Later  results  indicate  the  same  to  be  true 
for  all  configurations  utilizing  the  face  curtain  mode  of  egress.  For 
the  head  free  or  secondary  "D"  ring  configuration,  position  2,  the  maxi¬ 
mum  normal  stress  occurs  in  the  middle  and  upper  thoracic  region  due  to 
the  effects  of  head  movement  and  head  impact  on  the  anterior  chest  wall. 
Figure  21  shows  the  maximum  axial  force  and  bending  moment  distribution 
for  the  5th  percentile  aviator.  The  bending  moment  values  are  zero  for 
L5  and  L4  since  only  moments  tending  to  flex  the  column  were  considered. 
The  time  histories  of  axial  force,  bending  moment,  and  shearing  force 
for  LI,  the  level  of  maximum  normal  stress  for  the  5th  percentile  avia¬ 
tor,  are  given  in  Figure  22.  As  can  be  seen,  the  axial  force  level 
parallels  the  rise  of  the  acceleration  as  given  in  Figure  13.  Since  the 
axial  force  remains  relatively  constant  after  130  ms,  the  bending  moment 
dictates  where  the  maximum  normal  stress  occurs.  The  head  restraint 
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BENDING  MOMENT 


SHEAR  FORCE 


-1000  «- 


Figure  17.  Time  History  of  Axial  Force,  Shear  Force,  and  Bending  Mo 
ment  on  Superior  Surface  LI,  5th  Percentile,  Position  1, 
Unrestrained 


Figure  18.  Axial  Force  and  Bending  Moment  Distribution  at  170  ms, 
5th  Percentile,  Position  1,  Unrestrained 
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Figure  19.  Spinal  Alignment,  5th, 
Face  Curtain,  Shoulder 
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40th,  95th  Percentile,  Position  1 
Harness 


Figure  20.  Maximum  Normal  Stress,  5th,  40th,  95th  Percentile,  Posi 
tion  1,  Face  Curtain,  Shoulder  Harness 
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AXIAL  FORCE,  lb;  SHEAR,  lb 
BENDING  MOMENT,  in-lb 


Figure  21.  Maximum  Axial  Force  and  Bending  Moment  Distribution,  5th 
Percentile,  Position  1,  Face  Curtain,  Shoulder  Harness 


Figure  22.  Time  History  of  Axial  Force,  Shear  Force,  and  Bending  Mo¬ 
ment  of  Superior  Surface  of  LI,  5th  Percentile,  Position  1 
Face  Curtain,  Shoulder  Harness 
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Figure  23.  Face  Curtain  Force,  5th  Percentile,  Position  1,  Face 
Curtain,  Shoulder  Harness 
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force  is  shown  in  Figure  23.  Peak  values  of  this  force  are  in  the  same 
order  of  magnitude  as  the  results  of  Vulcan  and  King  (1970) . 

The  results  of  Configuration  3, shown  in  Figures  24  through  27, 
give  the  outputs  for  the  5th,  40th  and  95th  percentile  aviator  with 
face  curtain,  shoulder  harness  and  seat  back  constraints.  The  seat 
back  constraint  was  idealized  to  give  lateral  but  no  vertical  support. 
There  probably  exists  a  certain  amount  of  vertical  support  due  to 
friction,  as  is  the  case  with  the  shoulder  harness,  but  this  effect  was 
neglected.  Comparing  the  results  of  Figures  19  and  24  shows  that  the 
lateral  restraint  prevents  the  bowing  oif  the  lower  column.  Figure  25 
shows  a  75%  reduction  in  normal  stress  due  to  the  lateral  restraint. 

The  amplitudes  of  peak  normal  stress  have  also  generally  shifted  to  a 
level  lower  in  the  spinal  column.  The  5th  percentile  peak  occurs  at 
L2,  the  40th  at  T10,  and  the  95th  at  T10.  When  lateral  movement  is  cur¬ 
tailed,  the  greater  influence  of  axial  stress  on  the  normal  stress  is 
apparent.  The  axial  stress  comprises  approximately  30%  of  the  normal 
stress  with  lateral  support  whereas  with  no  lateral  support,  it  only 
represented  8%.  Hence,  one  can  conclude  that  lateral  support,  in  the 
form  of  a  properly  designed  back  pad  and  harness  system,  can  prevent 
excessive  anterior  compressive  stresses  with  its  resulting  anterior  lip 
fracture.  The  seat  back  itself  would  not  necessarily  provide  the  de¬ 
gree  of  support  required.  There  is  probably  a  substantial  degree  of 
movement  allowed  since  the  curvatures  of  the  spine  would  prevent  inti¬ 
mate  contact.  Additionally,  the  rib  cage  humps  posteriorly  and  pre¬ 
vents  spinous  processes  from  contacting  the  seat  back.  Since  a  lateral 
deflection  of  only  0.5  in.  is  sufficient  to  cause  large  bending  stres¬ 
ses,  a  pad  which  is  specifically  and  intimately  contoured  for  each  avi¬ 
ator  might  be  one  of  the  simpler  solutions  to  the  ejection  problem. 
Alternatively,  the  support  might  be  designed  to  place  all  spinal  regions 
in  extension  so  as  to  load  the  posterior  lamina  and  pedicles  and  thereby 
increase  the  effectively  loaded  area.  Figure  26  shows  the  maximum  axial 
force  and  axial  stress  levels.  Again,  it  can  be  seen  that  axial  stress 
levels  remain  fairly  constant.  The  increase  in  axial  force  levels  from 
the  5th  to  95th  percentile  aviators  is  due  to  the  increase  in  mass. 
However,  it  should  be  noted  that  the  cross-sectional  area  of  the  discs 
and  vertebrae  used  were  the  same  for  all  runs.  Since  disc  and  verte¬ 
bral  areas  increase  in  proportion  to  overall  body  size,  it  is  possible 
that  the  spines  of  each  aviator  would  be  subjected  to  the  same  axial 
stress  levels.  Thus,  the  initial  spinal  configuration,  through  its  in¬ 
fluence  on  the  location  of  maximum  bending  stress,  would  govern  the 
level  of  anterior  spinal  fracture.  The  time  histories  of  force,  bending 
moment,  and  shearing  force  at  L2  for  the  5th  percentile  aviator  are 
given  in  Figure  27. 

Figures  28  through  31  give  the  results  of  Configuration  4,  i.e., 
the  5th  percentile  with  shoulder  harness  and  seat  back  restraint,  but 
the  head  is  free.  The  displacement  of  the  Tl-6  segment  seen  in  Figure 
28  would  not  be  as  pronounced  if  the  column  had  been  discretized.  If 
this  had  been  done,  the  bending  and  axial  forces  would  have  been  dis¬ 
tributed  over  the  additional  discs  and  vertebrae.  However',  the 
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Figure  24.  Spinal  Alignment,  5th,  40th,  95th  Percentile,  Position  1 
Face  Curtain,  Shoulder  Harness,  Seat  Back 
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Figure  25.  Maximum  Normal  Stress,  5th,  40th,  95th  Percentile,  Posi 
tion  1,  Face  Curtain,  Shoulder  Harness,  Seat  Back’ 
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Figure  27.  Time  History  of  Axial  Force,  Shearing  Force  and  Bending 

Moment  on  Superior  Surface  of  L2,  5th  Percentile,  Position 
1,  Face  Curtain,  Shoulder  Harness,  Seat  Back 


LATERAL  DISPLACEMENT,  in 

Figure  28.  Spinal  Alignment,  5th  Percentile,  Position  2,  Shoulder 
Harness,  Seat -Back 


AXIAL  FORCE  ,  lb  ;  SHEARING  FORCE  ,  lb  and 


Figure  29.  Maximum  Normal  Stress,  5th  Percentile,  Position  2,  Shoulder 
Harness,  Seat  Back 


Figure  30.  Time  History  of  Axial  Force,  Shearing  Force,  and  Bending 
Moment  on  Superior  Surface  of  T7,  5th  Percentile,  Posi¬ 
tion  2,  Shoulder  Harness,  Seat  Back 
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AXIAL  FORCE,  lb 
BENDING  MOMENT,  in-lb 


'Figure  31.  Maximum  Axial  Force  and  Bending  Moment,  5th  Percentile, 
Position  2,  Shoulder  Harness,  Seat  Back 


LATERAL  DISPLACEMENT ,  in 

Figure  32.  Spinal  Alignment,  5th  Percentile,  McDonnell  Lumbar  Pad, 
Face  Curtain,  Shoulder  Harness,  Seat  Back 
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configuration  does  serve  to  show  that  the  maximum  normal  stress  is 
shifted  to  the  upper  thoracic  region  (Figure  29)  when  the  face  curtain 
is  not  used.  It  also  explains  the  observation  of  Latham  (1957)  that 
inadvertent  ejection  with  loose  harness  and  failure  to  utilize  the  face 
curtain  always  resulted  in  a  fracture  in  the  region  of  T6.  The  time 
histories  of  axial  force,  bending  moment,  and  shearing  force  at  T7  are 
given  in  Figure  30.  Figure  31  gives  the  maximum  axial  force  and  bend¬ 
ing  moment  distribution  for  the  spinal  column  with  the  head  free.  In 
each  figure  the  effect  of  the  freely  displacing  head  can  be  seen. 

One  final  configuration  was  studied.  Sanford  and  Kellet  (1966) 
reported  on  an  improved  lumbar  pad  to  be  used  by  the  U.S.  Navy  in  the 
MK-H5/H7F-1  seat.  The  lumbar  curvatures  reported  for  a  15th  percentile 
aviator  were  used  along  with  the  anthropomorphic  data  for  the  5th  per¬ 
centile  aviator  as  reported  by  Kaplan  (1972).  Figures  32  through  35 
give  the  results  for  the  face  curtain,  shoulder  harness  and  seat  re¬ 
straint  configuration.  The  main  effect  of  the  lumbar  pad  is  to  place 
the  column  in  hyperextension  as  can  be  seen  in  Figure  32.  This  allows 
the  posterior  vertebral  elements  to  carry  a  portion  of  the  axial  and 
bending  load  and  thus  greatly  decreases  the  anterior  stress  as  can  be 
seen  in  Figure  33.  However,  the  axial  force  levels  and  stresses  are 
increased.  Figure  34  shows  the  time  history  of  the  axial  force,  bend¬ 
ing  moment  and  shearing  stress  at  T7.  Figure  35  gives  the  axial  force 
and  bending  moment  distribution  at  150  ms. 

A  synopsis  of  the  results  for  the  various  runs  is  given  in  Table 
7.  Maximum  axial  forces  are  seen  to  occur  in  the  lower  lumbar  region. 
However,  since  disc  cross-sectional  areas  vary,  the  maximum  axial 
stresses  occur  in  the  lower  thoracic  and  upper  lumbar  regions.  Thus, 
initial  configurations  which  place  the  level  of  maximum  bending  stress 
in  the  thoracolumbar  region  should  be  avoided,  since  the  maximum  normal 
stress  then  indicates  a  strong  probability  of  anterior  lip  fracture 
there.  Yamada  (1970)  found  compressive  vertebral  failure  occurred  when 
the  body  of  the  vertebra  was  exposed  to  stresses  greater  than  908  to 
1249  lb/in2.  Values  in  Table  7  indicate  that  with  inadequate  lateral 
support  stress  levels  of  up  to  200  to  300%  of  the  failure  level  can 
occur.  If  the  head  is  free,  the  stress  levels  obtained  indicate  a  very 
strong  likelihood  of  fracture.  If  one  accepts  the  use  of  effective 
normal  stress  as  a  criterion,  the  necessity  of  utilizing  the  face  cur¬ 
tain  and  having  adequate  seat  back  support  is  very  clearly  indicated. 

The  fact  that  the  situation  is  not  quite  as  critical  as  indicated 
in  the  proceeding  paragraph  is  probably  due  to  two  sources: 

1.  One  of  the  basic  assumptions  of  the  Orne-Liu  model  is  the 
idealization  of  the  soft  tissues  as  a  rigid  body.  Thus,  of  the  two  as¬ 
signed  properties  of  the  soft  tissues,  inertia  and  elasticity,  only 
the  inertia  has  been  retained.  .  This  is  a  cautious  formulation  of  the 
problem  since  the  deformations  of  the  soft  tissues  would  alleviate  to  a 
certain  extent,  the  stresses  experienced  by  the  vertebrae.  This  ap¬ 
proach  was  adopted  to  make  the  analytic  problem  tractable  while 
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VERTEBRAL  LEVEL 


Figure  34.  Time  History  of  Axial  Force,  Shearing  Force,  and  Bending 

Moment  on  Superior  Surface  of  T7,  5th  Percentile,  Position  1, 
Face  Curtain,  Shoulder  Harness,  Seat  Back,  McDonnell  Lumbar 
Pad 
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Figure  35. 


Axial  Forces  and  Bending  Moments  at  T  =  150  ms,  5th  Per¬ 
centile,  Position  1,  Face  Curtain,  Shoulder  Harness,  Seat 
Back,  McDonnell  Lumbar  Pad 


TABLE  7 


TEST  RESULTS 


CONFIGURATION 

PERCENTILE 

AVIATOR 

MAXIMUM 
AXIAL  FORCE 
lb 

VERT. 

LEVEL 

TIME 

ms 

MAXIMUM 
AXIAL  STRESS 
psi 

VERT. 

LEVEL 

TIME 

ms 

MAXIMUM 

BENDING 

MQMENT+ 

VERT. 

LEVEL 

time 

ms 

MAXIMUM 
NORMAL  STRESS 
psi 

VERT. 

LEVEL 

TIME 

ms 

Unrestrained 

5 

457 

L5 

130 

301 

Til 

150 

1065* 

LI 

170 

7211 

LI 

170 

FC,  SH 

5 

610 

L5 

150 

260 

Til 

150 

561 

LI 

170 

3946 

LI 

170 

40 

698 

L5 

150 

282 

L2 

150 

653 

T12 

200 

4555 

T12 

200 

95 

878 

L5 

150 

402 

L2 

150 

477 

T9 

180 

4239 

T9 

180 

FC,  SH,  SB 

5 

746 

L5 

150 

346 

T9 

160 

91 

L2 

160 

832 

L2 

160 

40 

840 

L5 

160 

369 

L2 

160 

136 

T10 

170 

1286 

T10 

170 

95 

959 

LS 

140 

428 

L2 

140 

158 

T10 

160 

1463 

T10 

160 

Head  Free,  SH, 
SB 

5 

711 

L5 

140 

156 

T9 

130 

1416 

T7 

210 

15713 

T7 

210 

FC,  SH,  SB, 

MLP 

5 

859 

L5 

150 

349 

Til 

150 

98* 

T7 

150 

1365 

T7 

ISO 

FC  =  Face  Curtain 
SH  =  Shoulder  Harness 
SB  =  Seat  Back 
MLP  =  McDonnell  Lumbar  Pad 

t  in- lb 

•  Secondary  Maximum  chosen  since  LS  was  clamped. 


providing  an  upper  bound  estimate  of  the  force  levels  found  in  the  real 
case . 


2.  The  failure  loads  given  by  Yamada  (1970)  are  from  quasistatic 
tests.  It  is  well-known  that  at  the  high  strain  rate  encountered  in  pi¬ 
lot  ejection,  the  fracture  stress  is  considerably  higher  than  indicated 
by  these  values. 


CONCLUSIONS 


1.  Spinal  alignment  prior  to  ejection,  i.e.,  initial  configuration,  is 
■  a  major  factor  in  determining  the  location  and  magnitude  of  maximum 

vertebral  stresses  for  a  given  set  of  restraints. 

2.  A  properly  designed  back  pad  and  harness  system  could  prevent  ex¬ 
cessive  anterior  vertebral  stresses  when  the  face  curtain  is  util¬ 
ized. 

3.  The  secondary  "D"  ring  ejection  position  causes  excessive  normal 
stresses  on  the  anterior  edges  of  upper  and  middle  thoracic  verte¬ 
brae  . 

4.  The  face  curtain  and  shoulder  harness  restraints  serve  to  place  the 
maximum  normal  stress  in  regions  of  maximum  spinal  curvature. 

5.  The  locations  of  maximum  normal  stress  levels  can  be  predicted  from 
spinal  curvatures  as  determined  from  x-rays  under  static  conditions 
for  the  face  curtain  ejection  position. 

6.  The  occurrence  of  maximum  normal  stress  corresponds  to  the  occur¬ 
rence  of  maximum  bending  moments. 

7.  The  effect  of  seat  back  restraint  is  to  reduce  maximum  normal  stres¬ 
ses  and  shift  their  location  caudally. 

8.  A  strong  probability  of  anterior  lip  failure  exists  when  the  face 
curtain  is  not  utilized  or  if  inadequate  (posterior-anterior)  sup¬ 
port  is  given  to  the  spinal  column. 

9.  The  use  of  a  lumbar  pad  which  places  the  lumbar  spine  in  extension 
greatly  reduces  the  anterior  normal  stresses. 

10.  Maximum  axial  stresses  occur  in  the  lower  thoracic  and  upper  lumbar 
regions.  Secondary  maximums  occur  in  the  middle  and  upper  thoracic 
column . 

11.  Variations  in  the  cross-sectional  areas  of  the  discs  and  vertebrae 
result  in  a  state  of  relatively  constant  axial  stress  being  main¬ 
tained  in  the  spinal  column  during  ejection. 
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RECOMMENDATIONS 


1.  Any  ejection  seat  design  or  improvements  should  be  studied  in  a  math¬ 
ematical  model  such  as  used  here  to  evaluate  the  effectiveness  of  the 
design  or  improvement  in  reducing  the  intervertebral  joint  stresses. 

2.  An  experimental  program  is  urgently  needed  to  determine  the  stress 
levels  required  and  the  mechanisms  involved  in  vertebral  failure  un¬ 
der  combined  states  of  stress.  Proper  design  of  restraint  and  sup¬ 
port  systems  should  be  predicated  on  this  knowledge. 

3.  Studies  should  be  undertaken  to  provide  a  statistically  significant 
set  of  spinal  alignment  data  for  pilots  with  a  view  toward  develop¬ 
ing  a  percentile  spinal  configuration  population.  To  complement 
this  data,  anthropometric  measurements  needed  for  inertial  proper¬ 
ties  determination  should  be  gathered  concurrently,  i.e.,  height, 
weight,  sitting  height,  chest  depth,  chest  breadth,  and  hip  breadth. 

4.  Upon  completion  of  the  experimental  program,  a  spinal  support  sys¬ 
tem  should  be  designed  to'  minimize  stresses  during  ejection.  This 
should  be  accomplished  using  data  gathered  during  the  experimental 
phases  and  available  mathematical  modelling  techniques. 

5.  Any  improved  restraint  design  should  be  capable  of  being  installed 
on  existing  ejection  seats. 

6.  Complete  spinal  alignment  x-rays  should  be  taken  while  the  McDonnell 
lumbar  pad  and  any  other  similar  devices  are  in  place.  The  data 
should  then  be  subjected  to  analysis  by  currently  available  techni¬ 
ques  . 
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